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ABSTRACT 


Analyses of multlspectral measurements from satellites offer the 
potential to monitor and inventory crop production. The crop canopy is 
a dynamic entity influenced by many cultural and environmental factors. 

In order to quantify and understand these potential sources of variation 
in spectral measurements of crops , an experiment was conducted at the 
Williston, North Dakota, Agricultural Experiment Station in 1977. The 
effects of soil moisture, planting date, nitrogen fertilization, and 
cultivar on reflectance of spring wheat (Triticum aestivum L .) canopies 
were investigated. Spectral measurements were acquired on eight dates 
throughout the growing season with a radiometer (Exotech 100 A) in four 
wavelength hands (0.5-0. 6, 0.6-0. 7, 0.7-0. 8, 0.8-1. 1 pm) at 3.5 m above 
the plots. On each date that spectral data were acquired, concomitant 
measurements of crop maturity stage, leaf area index, biomass, plant 
height, percent soil cover, and soil moisture were made. 

In this experiment, planting date and available soil moisture were 
t’te primary agronomic factors which affected reflectance of spring wheat 
canopies from tillering to maturity. Comparions (R ) of treatments indi- 
cated that during the seedling and tillering stages planting date was 
associated with 36 and 85 percent of variation in red (0.6-0. 7 pm) and 
near Infrared (0. 8-1.1 pm) reflectances, respectively. As the wheat 
headed and matured, less of the variation in reflectance was associated 
with planting date and more with available soil moisture. By mid July 
soil moisture accounted for 73 and 69 percent of the variation in 
reflectance in red and near infrared bands, respectively. Differences in 
spectral reflectance among treatments were attributed to changes in leaf 
area index (LAI), biomass and percent soil cover. Cultivar and N fertiliza- 
tion rate were associated with very little of the variation in reflectance 
of these canopies. Agronomic practices which result in differences in LAI, 
biomass, and percent soil cover potentially can be monitored by remote 
sensing and may be useful in estimating crop production. 


Contribution from the Laboratory for Applications of Remote Sensing and 
Department of Agronomy, Purdue University. Journal Paper 8054, Purdue 
Agriculture Experiment Station, West Lafayette, IN. 47907. This study was 
supported by National Aeronautics and Space Administration Contracts 
NAS9-14970 and NAS9-15466. 
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INTRODUCTION 


The reflectance properties of single plant leaves have been identified 
and studied in laboratories for more than a decade. Relationships of physical 
and biological characteristics, such as leaf morphology, chlorophyll concen- 
tration, and leaf water content, to reflectance of leaves have been well 
established. Several comprehensive reviews of various aspects of this 
research have been published (2,4,5,7,12). 

Although knowing the reflectance characteristics of single leaves is 
basic to understanding the reflectivity of crop canopies, significant dif- 
ferences exist between spectra of single leaves and crop canopies. Canopies 
are more than simple collections of leaves. Complex interactions which are 
not factors when spectra of single leaves are measured must be considered in 
remote sensing of canopies grown under field conditions. Some important 
agronomic parameters influencing reflectance of crop canopies are leaf 
area index, percent soil cover, biomass, leaf geometry, leaf color and soil 
color (2,3,7). Because the crop canopies are dynamic entities, they are 
influenced by many management practices and environmental factors, including 
cultivar seeding rate, fertilization, soil moisture and disease (1,2,3,9,14). 

Crop identification and crop area estimation promise to be two of the 
major applications of remote sensing. In the recently completed Large Area 
Crop Inventory Experiment (LACIE), this technology was pushed to near opera- 
tional use for wheat (10). Remote sensing also offers considerable potential 
for acquiring information about crop conditions and yields (2,14). If multi- 
spectral remote sensing is to be used successfully to Identify and inventory 
crops, it is important to quantify and understnad the sources of variation in 
spectral measurements of crops. Some of the variation may be associated with 
important agronomic factors, which it may be desirable to monitor or inventory 
(e.g. , dryland vs. irrigated wheat). On the other hand, it is also important 
to know the magnitude of tariation associated with a factor, such as cultivar, 
which one probably would not be interested in monitoring. 

The objectives of this experiment were (1) to determine the relationship 
of canopy characteristics to the reflectance of spring wheat and (2) to 
determine the effect of several agronomic treatments, which represented 
the farming practices of the northern U.S. Great Plains, on the multispectral 
reflectance of spring wheat during different development stages. 


MATERIALS AI^D METHODS 


Data x^^ere collected at the North Dakota State Agricultural Experiment 
Station at V7illiston (43,32°N, 103,42°W) during the summer of 1977. The 
station is located in the gently rolling uplands above the Missouri River 
Valley and is representative of dryland farms of the region. Because 
of limited precipitation (36 cm per year) the majority of the land is planted 
to crops every other year and is left fallow in intermediate years to 
accumulate subsoil moisture. 


The spring wheat ( Triticum aestivum L.) experiment was a split plot 
design. Within each available soil moisture condition (whole plot), 
there were two blocks replications of a factorial experiment with cultivar, 
nitrogen fertilization, and planting date as treatments: 


Soil Moisture: 
Planting Date: 
Cultivar: 
Nitrogen Fertilization: 


Low (wheat crop in 1976) 
High (lallow in 1976) 

Early (9 May 1977) 

Late (23 May 1977) 

Semi-dwarf (Olaf) 

S tandard (Waldron) 

None 

44 kg N/ha 


The plots were 3.5 m wide and 15.3 m long with 18 cm wide, north-south rows. 
The soil was a Williams loam fine-loamy, mixed typic argiborolls which has 
a dark brown (10 YR 3/2) surface color when moist and very light (10 YR 4/3) 
color when dry. Although moisture in the top 20 cm of soil of the whole 
plots were similar at planting, the profile of fallow soil contained 20 per- 
cent (3 cm) more water in the 20 to 60 cm zone than the profile of soil on 
which wheat was grown in the previous year. 

Spectral measurements were made with an Exotech 100 radiometer in four 
wavelength bands, 0.5-0. 6, 0.6-0. 7, 0. 7-0.8, and 0. 8-1.1 ym, corresponding 
to the four spectral bands of the Landsat MSS. Measurements in all bands were 
taken concurrently and recorded by a printing data logger. Duplicate obser- 
vations were acquired over each plot and were averaged for these analyses. 

A boom mounted on a van supported the radiometer vertically at 3.5 m above the 
canopies and 3.5 m away from the van. At this elevation with 15° field of 
view the sensor viewed a 0-9 m diameter ground area. The reflectances of all 
plots were measured on eight dates during the growing season when the sun 
angle was greater than 45° above the horizon. On selected dates reflectances 
were measured four to eight time per day at approximately hourly intervals. 
Only data acquired during a two hour interval centered on solar noon were 
used in these analyses. 


A 1.2 m square ranel painted Vxth highly reflecting barium sulfate 
was used as a reference surface for determination of the bidirectional 
reflectance factor (11) . This reflectance standard provided a field 
calibration reference with stable, known reflectance properties. A 
dark level response of the instrument was also obtained by holding an 
opaque, light-ti,'>ht apparatus against the instrument's optical ports 
to measure the Internal system noise or deviation from zero. The 
response of the reference panel was measured approximately every 15 
minutes during the data collection period and the dark level every 30 
minutes. These values were then used in the following equation to calibrate 
readings taken over the plots: 

{xf • 


Where, BRF (A) * bidirectional reflectance factor (%) at a specific 
wavelength Interval (X) 

Ds(A) = response of Instrument to scene (plot), 

ds(A) = dark level response of instrument taken closest to 
time to scene , 

Dr (A) = response of instrument to painted barium sulfate 

reference standard, 

dr(A) = dark level response of instrument taken closest in 
time to reference standard measurement, and 

Rr(A) « reflectance of painted barium sulfate reference 

standard (measurement in laboratory by comparison 
with pressed barium sulfate) . 

Spectral data were acquired over the north end of these plots and 
agronimic data were collected from the southern two-thirds of each plot. 

This kept the plants viewed by the sensor untrampled and intact throughout 
the growing season. Measurements of crop maturity stage (8), percent soil 
cover, plant height, leaf area index, biomass, soil moisture, and 
comments on crop condition were acquired each day that reflectance data 
were acquired. Crop biomass in 2 each plot was estimated by harvesting plants 
in 1.0 m length of row (0.178 m ). Each sample was placed in a plastic bag; 
weighed (fresh biomass); separated into stems, heads, green, yellow, and brown 
leaf blades (leaf sheaths were included with stems); dried at 60°C and 
reweighed. The leaf area of a random subsample of green leaf blades from 
each plot was measured (Lambda Instrument Co., Model LI-3000) and the 
leaf area: leaf dry weight ratio was estimated. Leaf area index was 

calculated using this leaf area: leaf weight ratio and total dry weight 

of green leaf blades from plants in 1.0 m length of row. Grain yields 
were estimated from a 7.4 m area harvested from the center of each plot with 


a small self-propelled combine. Vertical color photographs taken from 6 ra 
were used to estimate percent soil cover in each plot. 

Dally meteorological data useful for describing the growing season 
were acquired at a National Weather Service cooperative station located on 
the experiment station. On each day that spectral data were collected 
additional meteorological measurements including air temperature, barometric 
pressure, relative humidity, wind speed and direction, and total irradir.nce 
were recorded continuously on strip charts. These data were used primarily 
to document atmospheric conditions on days when spectral data were acquired. 

The reflectance data were analyzed as band means and as transformations. 
The reflectance data were transformed using a principal component analysis 
into "greenness" and "brightness" variables as described by Kauth and Thomas 
(6) for Landsat MSS data and modified for spectrometer data. Correlation 
analyses were used to quantify the relationship between spectral values and 
agronomic characteristics. Analysis of variance was used to determine which 
of the experimental treatments accounted for the variability in spectral 
responses . 
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RESULTS AND DISCUSSION 


Relation of Canopy Characteristics to Reflectance 


The amount of vegetation present Is one of the principal factors affect- 
ing the reflectance of crop canopies. Figure 1 Illustrates the relationship 
between an agronomically Important canopy characteristic, leaf area Index 
(LAI) , and reflectance In selected wavelength bands . This figure Includes 
data from all treatments when green leaves were present (from seedling through 
heading). A portion of the scatter in the data is associated with various 
agronomic treatments, as well as measurement errors in the independent and de- 
pendent variables. 

As LAI Increased, red (0.6-0. 7 pm) reflectance decreased and near 
Infrared (0. 8-1.1 pm) reflectance Increased (Figure 1). These relationship 
of LAI and reflectance are slightly nonlinear, particularly In the red band* 
Studies of other canopies have Indicated asymptotic responses of reflectance 
at leaf area Indices greater than 3 or 4 (3.14). 

Correlations of reflectance measurements in Landsat radiometer bands 
with several agronomic characteristics of spring wheat canopies are shown 
In Table 1. Fresh biomass, dry biomass, and plant water content correlated 
most highly with reflectance in the visible wavelengths and with the greenness 
transformation. Leaf area Index and percent soil cover correlated highly 
with red and near Infrared reflectance and the greenness trans format:* on. 

Previous research has indicated that the amount of photosynthetically active 
(green) vegetation was highly correlated with reflectance of crop canopies 
(1,3,9). The decreased correlations of canopy variables and reflectance as 
the wheat began to ripen and senesce substantiated these observations (Table 1). 

Plant growth and development were significantly Influenced by the level 
of available soil moisture (Table 2) and these changes in canopy characteristics 
were manifested in their reflect£Uices (Figures 2 and 3). For example, canopies 
of wheat planted early on fallow soil (high soil moisture) had sign!. Icantly 
greater biomass, LAI, percent soil cover, infrared reflectance and lower red 
reflectance than canopies of other treatment combinations (Table 2). 

Precipitation on the day prior to the day when spectral data were acquired 
tended to decrease both red and Infrared reflectances. Approximately 10 
to 12 ram of rain, which darkened the soil, fell prior to the acquisition of 
spectral measurements on 23 June and 7 July and probably contributed to the 
abrupt decreases in reflectance on these dates. These decreases are most 
evident for those canopies of late planted wheat (Figure 2 and 3) which had 
the lowest soil cover percentages (Table 2). Thompson (13) noted similar 
decreases in scene radiances following precipitation and was able to de- 
lineate precipitation patterns and drought severity using Landsat MSS Imagery. 
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Leaf Area Index 

Figure 1. Relationships of red (0.6-0. 7 ym) and near infrared (0. 8-1.1 
ym) reflectance and leaf area index of spring wheat canopies 
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Table 1. Correlations of bidirectional rofJcc taiif c l.n'tn wIfU n 

agronomic charactarlstlcf} of spring wheat i if. fi » fh't f 
time intervals. 


Wavelength 
Band, urn 

Plant 

Height 

Percent 

Soil 

Cover 

Leaf 

Area 

Index 

Fresh 
B1 oina.sr 

ih V 

!* 



Seedling to 

Harvest 

Maturity 



0.5-0. 6 

-0.40 

-0.77 

-0.74 

“0.74 

“0.54 

- (i. 76 

0.6-0. 7 

-0.23 

-0.69 

-0.79 

"0.67 

“0.36 

0 75 

0. 7-0.8 

0.20 

0.14 

0.52 

0.13 


Lh 32 

0. 8-1.1 

0.02 

0.45 

0.81 

0.s5 

0.U3 

0.64 

Greenness 

0.17 

0.67 

0.88 

0.65 

0.25 

.81 

Sample size 

256 

224 

165 

224 

22a 

2 74 
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Seedling 

to Soft Dough 



0.5-0. 6 

-0.63 

-0.83 

-0.79 

-0.80 


-0.84 

0.6-0. 7 

-0.55 

“0.80 

-0.81 

-0.74 

^ 0.57 

-0.80 

0. 7-0.8 

0.03 

0.30 

0.52 

0.23 

0.03 

0.32 

0. 8-1.1 

0.37 

0.66 

0.81 

0.61 

0-3‘J 

0 70 

Greenness 

0.51 

0.80 

0.88 

0.75 

0.53 

0.82 

Sample size ■ 

153 

153 

153 

353 

3 53 

153 




Seedling 

to Flowering 



0. 5-0.6 

-0.76 

“0.89 

-0.86 

“0.85 

“0.76 

- U .86 

0.6-0. 7 

-0.73 

-0.88 

-0.86 

“0.82 

-0. 73 

-0.84 

0.7-0. 8 

0.37 

0.51 

0.56 

0.46 

0.34 

0 49 

0. 8-1.1 

0.65 

0.80 

0.83 

0.76 

0.64 

0.79 

Greenness 

0.75 

0.91 

0.91 

0.85 

0.74 

0.88 

Sample s ize » 

101 

301 

101 

103 

iOj 

101 


.Plant water content “ 


Fresh Biomass - Dry Blouas.*; 


Mean agronomic characteristics of spring wheat canopies as influenced 
available soil moisture levels and two planting dates. 
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Figure 2. Temporal changes in red (0.6-0. 7 pm) bidirectional 

reflectance factor of Spring wheat canopies at two levels 
of available soil moisture and two planting dates. The 
standard errors of the mean on each date are smaller than 
the symbols used on these graphs. Occurrence and amount of 
precipitation are also indicated. 


10 


Precipitation (MM/day) Bidirectionai Reflectance Factor (%) 


Eariy Planting Date 

(May 9,1977) 


Late Planting Date 

(May 23,1977) 


Soil Moisture 

High (Fallow In 1976) 
• Low (Wheat in 1976) 
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Figure 3. Temporal changes in near infrared (0. 8-1.1 pm) bidirectional 
reflectance factor of spring wheat canopies at two levels of 
available soil moisture and two planting dates. The standard 
errors of the means on each date are smaller than the symbols 
used on these graphs . Occurrence and amount of precipitation 
are also indicated. 



Relation of Agronomic Treatments to Reflectance 


The high correlations of red (0.6-0. 7 pm) and near infrared (0. 8-1.1 
pm) reflectance to green leaf area Index and percent soil cover (Table 1) 
potentially can provide the basis to inventory indirectly certain cultural 
and manageoent practices which could be useful In assessing crop production. 

For example If the ratios of near Infrared to red (IR/red) reflectances of 
these wheat canopies at heading are plotted against grain yields, two distinct 
groupings of data are observed (Figure A ) . Those wheat canopies with high 
IR/red reflectance also had high grain yields. Wheat grown on fallow soil 
(high soil moisture) produced 1748 and 1488 kg/ha of grain fur early and late 
planting dates, respectively, compared to 648 and 724 kg/ha of grain for 
early and later-planted wheat, respev':tlvely , with low soil moisture. 

Although no statistical test for the main effect of soli moisture 
level (whole plot in a nested design' was available, significant interactions 
of other factors with soil moisture provided Indications of its significance. 

A method of ascertaining the overall importance of the i, '•rimental treatments 
was to consider the percentage of total variability from .i,.alysis of variance 
(ANOVA) accounted for by each treatment in several wavelength bands and trans- 
formation (Tables 3, 4 and 5). Tnis approach permitted an evaluation of the 
effect of available soil moisture on reflectance of these canopies. 

On 1 June, planting date was the primary factor affecting spectral re- 
sponse (Table 3, 4 and 5). With a 14-day difference between the early and late 
planting dates there were significant differences In the agronomic, as well 
as spectral characteristics of the wheat canopies (Table 2). The wheat 
planted at the later planting date 23 May was at the seedling stage, while 
the early planted wheat (9 May) was beginning to tiller. The differences in 
spectral response may be partially attributed to the significant differences 
in percent soil cover, leaf area Index, and biomass, but since all wheat canopies 
covered less than 10 percent of the soil, differences in soil surface roughnw*'?s 
are probably the major contributors to the differences in spectral response. 
Spectral data were acquired after the soil of the early planted wheat had been 
smoothed and crusted by several rains, while the soil of the late planted 
wheat was still rough from planting. 

By early July differences in vegetative growth and spectral response 
due to planting date were reduced as the late planted wheat was nearly 
the same size as the early planted wheat (Table 2). At this time, heading 
to flowering stages of maturity, soil moisture availability was the major factor 
affecting the growth and development and spectral response of the wheat 
canopies . 

Although there were significant differences in July due to planting 
date, cultivar and nitrogen fertilization, they were small compared to those 
associated with available soil moisture. The primary difference in these 
two cultivars was the 'Waldron' reached its full height 2 to 3 days before 
'Olaf. Both cultivars developed similar leaf area and biomass. Ho'.’ever, 
we would expect that cultivars of wheat or any other crop which differ greatly 


Table 3. Percent oE variation (R ) in reflectance of 0.6-0. 7 um «ravelength band associated with 
soil tnoisturu (M) , cultivar(C) , nitrogen fertilization(.'0 • planting date V.P) and tlicir 
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Model Includes variation due to treatments. Total variation is due to blocks, 
treatments, and experimental error. Percentages less than 1.0 are omitted for 
clarity but are included in model. 
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Table 5. Percent of variation (R ) In the greenness transfoncatlon of reflectance associated with 
soil moisture (M) , cultlvar (C) , nitrogen fertilization (N) , planting date (F) and 
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In maturity and/or growth characteristics would also differ spectrally as 
shown by Learner et al. (9) • The addition of nitrogen fertilizer had only 
minor effects on the growth of wheat and accounted for very little of the 
variation in reflectance of these wheat canopies in 1977 when precipitation 
for June and July was below normal. 

On July 20, moisture level was the most Important factor accounting 
for differences in spectral response (Tables 3,4, and 5). Effects due to 
cultivar and fertilization were not significant except in the greenness t 
transformation (Table 5). 

On August 5, as the canopies were ripening, soil moisture level and 
planting date continued to be the primary factors influencing the spectral 
response especially in the 0.6-0. 7 pm band (Table 3) and greenness transformation 
(Table 5). Since the canopies had similar maturity stages, the difference 
in spectral response were attributed to the differences in percent soil cover 
and biomass (Table 2) . 

In summary, planting date and soil moisture status ware the primary 
agronomic factors affecting the reflectance of spring wheat canopies from 
tillering to maturity in a semi- arid environment. Changes in canopy reflectance 
due to planting dace and soil moisture were attributed to differences in maturity 
stage, leaf area index, biomass, and percent soil cover. Relatively high 
correlations between these canopy characteristics and reflectance were found. 
Cultivar and N fertilization had little effect on either the agronomic or 
spectral reflectance characteristics of these spring wheat canopies. Agronomic 
practices which result in differences in LAI, biomass, and/or percent soil 
cover can be monitored by remote sensing and may be useful in estimating 
crop production for large geographic areas. 
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